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Abstract 

Semiclassical oscillation of the electron through the nucleus of the H atom 
yields both the exact energy and the correct orbital angular momentum for 
/ = quantum states. Similarly, electron oscillation through the nuclei of 
H2 accounts for a stable molecule ion with energy close to the quantum 
mechanical solution. The small discrepancy arises from the neglect of the 
electron's wave nature. 
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I. INTRODUCTION 



Two of the reasons why the old quantum theory of Bohr and Sommerfeld 
was abandoned in the mid 1920s were the theory's failure to give the correct 
multiplet structure of the hydrogen atom and the stability of the hydrogen 
molecule ion, The old vector model of angular momentumP gave, for 

a given principal quantum number n, sublevels with angular quantum num- 
bers I = 1, 2, nP Spectroscopic evidence, however, showed multiplicities 
of spectral line splitting in a magnetic field according to / = 0, 1, ...,n — 1. 
Max Jammer, in his review,'^ notes that "the old quantum theory could never 
resolve this inconsistency." 

A treatment of the hydrogen molecule ion with Sommerfeld's quantization 
conditions had been Wolfgang Pauli's doctoral thesis of 1922.'^Pauli found its 
molecular binding energy to be positive (non-binding) — contrary to (later) 
experimental findings. Martin Gutzwillei'^ thinks that "the solution of this 
problem can be rated, with only slight exaggeration, as the most important 
in quantum mechanics, because if an energy level with a [more] negative 
value [than of a free hydrogen atom] can be found, then the chemical bond 
between two protons by a single electron has been explained." 

Both dilemmas of the old quantum theory can be resolved, though, with 
a single extension: an oscillation of the electron through the nucleus (nuclei) 
of the atom (molecule) . In essence this solution is already formally included 
in Sommerfeld's theory of the hydrogen atonPbut was explicitly omitted by 
Sommerfeld and his school as being unphysical. The case in point, obtained 
with Sommerfeld's quantization conditions for radial and angular motion, 
is a quantum state with zero angular action, characterized by an angular 
quantum number / = 0. What is its orbit? 

The geometry of an nl Sommerfeld ellipse is given by its semimajor axis, 
dni = {j'B/Z)n^, and semiminor axis, hni = {rB/Z)n^^/l{l + 1)^ Here tb = 
/i^/47r^me^ is the Bohr radius in terms of fundamental constants, serving as 
an atomic distance unit, and Ze is the nuclear charge. An (n, 0) orbit is thus a 
line ellipse with its nuclear focus at one end and its empty focus at the other. 
This case was regarded as unphysical because of the electron's collision with 
the nucleus — an uncritical adaptation from celestrial mechanics. A closer 
inspection confirms that a line ellipse with terminal nuclear focus is indeed 
unphysical — but for quite a different reason! 
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Fig. 1. Partial trajectory of an extranuclear orbit XX' and of 
a penetrating orbit PP' through nucleus A^. The dotted line S 
shows the major symmetry axis of XX'. 

Leaving quantization conditions momentarily aside, what would happen 
if we continuously decrease an angular quantum number A while keeping the 
principal quantum number n constant? We then would get more and more 
slender ellipses with the same length of major axis, 2a„. By basic electric 
theory, the nuclear Coulomb potential outside a finite-size nucleus is 
given by the point potential as if all nuclear charge, +Ze, was concentrated 
at the center of the nucleus. This holds as long as the electron orbit stays 
outside the nucleus. Two borderline cases are illustrated in Fig. 1. For a very 
small A value, say < Ax <^ 1, we obtain a very slender elliptical orbit with 
partial trajectory XX' about the nucleus. Further decrease of A to Ap < Ax 
causes an intrusion of the electron into the finite nucleus of radius (see 
trajectory PN in Fig. 1). Once inside, at a distance r < from the center, 
then, by Gauss's law, only a fraction of the nuclear charge, Z'{r)e < Ze, 
acts on the electron via centripetal force.l^ Accordingly, the electron's exit 
trajectory NP' is no longer symmetric to its approach trajectory PN with 
respect to the major axis 5* of the (partial) ellipse XX'. In the extreme 
case of a head-on penetration of the nucleus, A = 0, there is no centripetal 
force at all! The electron will then, with almost constant speed, traverse the 
nucleus, continue, with decreasing speed, to the opposite turning point of its 
line orbit and revert its motion periodically. We want to call the electron's 
straight-line oscillation in the Coulomb potential of a finite-size nucleus a 
"Coulomb oscillator. " 



Fig. 2. Line orbit of a Coulomb oscillator with nucleus at origin 
and turning points at ±r. 



II. HYDROGEN ATOM 

For the formal treatment of the non-relativistic Coulomb oscillator we 
designate the z axis along the line orbit, with turning points at z = iF and 
nuclear position at 2; = (see Fig. 2). The electron's total energy E at 
position z must equal the potential energy at a turning point, 

^ 1 2 Ze^ 
This gives the electron's speed along the z axis. 
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Fig. 3. Axial electron speed v vs. position z of an electron in 
Coulomb oscillation for the ground state, n = 1, of a hydrogen 
atom H {Z = 1, solid curve) and a helium ion He~^ {Z = 2, 
dashed curve) . The area under one wing of each curve represents 
the radial action, A — Ih. 



Figure 3 displays its dependence on the axial position as a two-wing curve 
cusped at the nucleus. Atomic units (a.u.) are used, that is, the Bohr radius 
vb and the "Bohr speed" vb = 27re^//i = ac — the electron speed in the 
ground-state Bohr orbit of the H atom — with fine-structure constant a ~ 
1/137 and speed of light c. The curve's wing along the positive z axis, 
1^1 = r, gives the radial speed, Vr{r) — \v{z)\, necessary for Sommerfeld's 
radial quantization condition, 

^ Prdr = m ^ Vr{r)dr = n^h. (3) 

Here p.,. is the radial momentum, n.f. = 1, 2, ... is the radial quantum number, 
and h is Planck's quantum of action. Integration is over one period of the 



5 



radial motion, z = +T — > — > +r. Graphically, the radial quantization is 
illustrated in Fig. 3 by the area under one wing of the speed curve. For a 
line orbit, / = , the radial quantum number equals the principal quantum 
number, n = rir + I = rir- 

In order to express the radial quantization in terms of axial motion we 
employ a "fold-out factor," = 2, to compensate for the doubling of inte- 
gration range in the extension from the radial one-wing speed curve to the 
axial two-wing curve. The quantization is thus restated, 

1 f TTl f 

— <p Pzdz = — <p v{z)dz = n^h, (3') 
9 J 9 J 

with axial quantum number Uz — rir — n and integration over the axial 
double-wing range, z = +r — > — F +F. By symmetry we can restrict 
the axial action integral to one quarter of the oscillation, say z — +r ^ 0, 



Pzdz = — / v{z)dz = — -y2Zm / \ dz = nh. (4) 

9 Jr 9 Jr \ z r 



Here the electron's motion in the negative z direction is accounted for by 
the negative sign. Although the electron's speed through a point nucleus 
diverges, ^^(O) = oo, the action integral, Eq. (4), stays finite and determines 
the quantized amplitude F„ of the Coulomb oscillator. Graphically the am- 
plitude F„ must be such that it stretches the speed curve horizontally to 
the extent that the area under one wing, = nh, represents the quantized 
action. The analytic solution, derived in Appendix A, is 

r„ = 2'fn\ (5) 
Inserting Eq. (5) into Eq. (1) yields the quantized energy. 



^2g2 ^2 

En = ^ = ~r^^2/' (6) 



in terms of the Rydberg energy unit, Ry = 2Ti'^me^ /h? — 13.6 eV , and in 
agreement with the energy of the nth Bohr orbit. 

Note that Eq. (5) gives the amplitude of the nth Coulomb oscillator 
as twice the radius of the nth Bohr orbit or of the semimajor axis of an 
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nl Sommerfeld ellipse, r„ = = {rB/Z)n'^. For comparison, the time- 
average radial distance of a Kepler orbilP of major and minor semiaxes a 
and b, respectively, is {r)t = (3a^ — 6^)/(2a). A line ellipse (6 = 0) has then 
{f)t = fa- Applied to an nl Sommerfeld orbit,'^ its average size, {rni)t = 
{rB/Z)[3n'^ — l{l + l)]/2, is in agreement with the corresponding quantity 
from quantum mechanics,'^ (vni) = J ip*r'ipd?r. Thus the time-average radial 
distance of a Coulomb oscillator is (r„o)t = \{rB/Z)in?. For the ground state 
of the hydrogen atom, n = 1, this gives (rio)t = fr^, as is well-known from 
quantum mechanics."^ 

The concept of the electron's semiclassical Coulomb oscillation is consis- 
tent with the Fermi-contact term of hyperfine interaction for / = states, 
which arises from the presence of the electron inside the nucleus. This is 
familiar from quantum mechanic^ and can be interpreted semiclassically as 
a local- field effect .l^^l 

To be sure, the extension of Sommerfeld's theory by the Coulomb oscil- 
lator resolves the discrepancy of the old quantum theory with spectroscopy, 
mentioned above, only at the low end of angular quantum numbers, / = 0. 
The resolution at the high end — repeal of the circular Bohr orbit, / 7^ n — 
involves an analysis in terms of space quantization which is beyond the scope 
of the present study. 








-i-c 
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Fig. 4. Axial Coulomb oscillation in an H2 molecule ion between 
axial turning points ±C and through nuclei at ±c; perpendicu- 
lar oscillation between lateral turning points ±5 and through 
midpoint 0. 
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III. HYDROGEN MOLECULE ION 



A. Formalism 

A hydrogen molecule ion, H2, consists of two proton nuclei and one 
electron. We assume, in adiabatic approximation, the protons located at fixed 
positions z = ±c on the molecular axis (see Fig. 4). The term "Coulomb 
oscillator" denominates again the motion of a point electron, now along either 
the line through the protons {z axis) and with turning points ±C, or along 
the perpendicular line through the midpoint (y axis) with turning points ±B. 
The molecule ion's total energy E at any position on the axis, — C ^ ^ +C, 
must equal the potential energy at the turning point C, 



E = -mv^ - - r H = — \ . (7) 

2 \z + c\ \z-c\ 2c C + c C-c 2c ^ ^ 

For electron positions beyond the protons, z > c, this gives an electron speed 



, 2e / ^ C 



For positions between the protons, < 2; < c, the corresponding speed is 



, 2e / c C 



The speed expressions will be used in the action integral, 
1 /■ , 2m 





with outer contribution 



2m 

■^out — ~ 



and inner contribution 



/2m f^'^ 
p^dz = y v{z)dz = = Aout + Ain, (9) 



j Vout{z)dz + J Vout{z)dz 



(10a) 



An = I Vin{z)dz. (10b) 
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Here is a fold-out factor to be specified below. 

If the conditions are such that the electron swings along the z axis through 
the midpoint 0, that is, f j„(0) > 0, then there exists also an oscillation along 
the y axis with lateral speed m, having the same total energy, 

E = -mu^ H = H . (11) 

2 v^2^ 2c y/WT^ 2c ^ ^ 

Equal energy at the axial and lateral turning points, E{C) = E{B), Eqs. (7) 
and (11), determines the latters' geometric dependence. 



c^ 



5= + — -3c2. (12) 



Solving Eq. (11) for the lateral speed. 



_ 2e / 1 1 

^7^\l V7T^ ~ 7WT^' ^^^^ 

provides the integrand of the action integral over a lateral oscillation, 



2m 



-B 



B 



Subtraction of the protons' mutual repulsion from the total energy E of 
axial or lateral motion, Eqs. (7) and (11), gives the electronic energy. 



e2 



Eei = E- -, (15) 

in its dependence on the proton separation, R = 2c. This brackets the 
molecular problem with known atomic results, Eq. (6), in the limits of i? = cxo 
(free H atom, Z = 1) and R = (free He~^ ion, Z = 2). For those cases, as 
well as any proton-proton distance R between, we keep the action constant, 

A = + Ay = nh. (16) 

Equation (16) is the Einstein quantization conditiorP^' — a generalization of 
Sommerfeld's quantization over separable variables — where the quantum sum 
equals the sum of action integrals over topologically independent paths in 
phase space 
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Fig. 5. Axial speed v vs. position z in the ground state of an 
H2 molecule ion (bold curve) and, for comparison, of a free H 
atom (thin curve) . Circles indicate the axial the positions of the 
nuclei, here with a large separation, i? = 6 tb- 

Here we treat the molecule ion only in its ground state, n = 1. Analytic 
solutions of the action integrals, Eqs. (lOab) and (14), are complicated due 
to elliptic functions. We therefore integrate numerically and visualize the 
integrals by the area under the corresponding speed curves. The bold curve 
in Fig. 5 shows the axial electron speed v{z) for a far proton separation, 
R = 6 r^. The electron, in its semiclassical motion, then oscillates only about 
(and through) the right proton. The area under the speed curve, Eq. (Sab), 
proportionally represents the ground state's unity of action, m § v{z)dz = (ph, 
with a fold-out factor = 2 in analogy to the free-atom case, Eq. (4). The 
thin curve shows, for comparison, the axial electron speed in a free H atom — 
familiar from Fig. 3 — centered at the same proton position, +c. The pull 
from the left proton (at — c) on the oscillating electron can be seen by the 
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distortion of the speed function v{z) and the redistribution of the area under 
the curve. 



10- 

viz) ■ 

"u(y) : 




Fig. 6. Axial speed v vs. position z in the ground state of an H2 
molecule ion (M-shaped bold curve) and, for comparison, of a free 
H atom (A-shaped thin curve, centered at the right nucleus). The 
n -shaped dotted curve, centered at 0, shows the lateral speed u 
vs. the perpendicular position y in the molecule. Circles indicate 
the axial positions of the nuclei, here with a small separation, 
i? = 2 rs. 

When, with closer proton separation R, as in Fig. 6, the electron swings 
past the midpoint 0, then the single-cusp speed curve v{z) from Fig. 5 — 
akin in shape to letter A — becomes doublc-cusped (akin to letter M) and 
symmetric with respect to the bisector (y axis). Now there is also a lateral 
oscillation with speed u{y) having the same total energy E. The equality of 
E in both cases can be seen in Fig. 6 by the equality of axial and lateral speed 
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at the midpoint, v(0) — u{0) — a position where the electron experiences the 
same potential in either case. For convenience the lateral speed u{y), though 
perpendicular to the proton axis, is displayed in Fig. 6 together with the axial 
speed v{z). The lateral speed curve u{y), drawn dotted, is readily recognized 
by its dome shape (n ). At the bifurcation value of the proton separation, 
R, where the electron starts swinging though the midpoint 0, the axial speed 
curve v{z) changes from its one-centered A shape to a two-centered M shape. 
The area under the axial speed curve then abruptly doubles, M(i? — 6) ^ 
2A{R + S), upon a very small change in proton separation, 6 <^ R. In order 
to keep the action integral continuous at R, the sudden area doubling is 
compensated by a corresponding doubling of the fold-out factor from (f) — 2 
for R > R to (j) = 4 for R < R. The same fold-out factor, = 4, must be 
used for the lateral action integral Ay, Eq. (14), as will become clear shortly. 

The unity of action, A = Ih, is visualized again in Fig. 6. To this end we 
compare the right half of the bold M-shape curve of axial electron speed in 
the molecule with the thin curve A (if) of the axial speed in a free H atom 
positioned at the right nucleus, +c. Due to attraction from the left nucleus, 
the right wing of the M curve is smaller than that of A{H) by the area of 
lobe L. On the other hand, the left flank of the A[H) curve that extends over 
the negative z axis is smaller than the left wing of the lateral speed curve fl 
by the area of shce S. The area under the free-atom curve is then 

A{H) ^^M + L+^n-S. (17a) 

The areas of lobe and shce are comparable, 

L^S. (17b) 

When the tiny notch to the right of the saddle point of M is taken into 
account, then the approximations (17ab) become equations and combine to 

M + n = 2A(i/). (18) 

The area under both the axial and lateral speed curves is thus four times 
the area under one wing of the free-atom curve, M+n = 4 x iA(if). Since 
the latter represents one quantum of action, /i, the combined area M-|-n 
visualizes its double fold-out, (f) — 4. 

With very close proximity of the nuclei, i? — > 0, the crests of the M curve 
start merging while its saddle point, v{0), keeps rising. In the R = limit of 
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fusing nuclei the axial electron speed becomes that of a free He^ ion, M— >■ 
K{He^), familiar from Fig. 3. Concurrently, the lateral speed curve fl rises 
at its peak, m(0) = f (0), and narrows at its base until it, too, turns into the 
speed curve of the free He^ ion, fl A{He~^). In the R = limit the three 
curves merge, M(0) = n(0) = A(iJe+). 

The results of the Coulomb-oscillator approach will be compared with 
another semiclassical calculation of H2, by Strand and Reinhardt.'^ These 
authors, like Pauli,'^ separate the equation of motion in spheroidal coordi- 
nates, ^ = (r+ + r_)/2c , 7] = (r_|_ — r_)/2c and ip by virtue of the constants 
of the motion: total energy E, angular momentum M about the z axis, and 
a component of the bifocal Runge-Lenz vector, Here r+ (r_) is the 

distance of the electron from the nucleus at +c (— c). Strand and Reinhardt 
(SR) solve the ensuing one-dimensional differential equations with classical 
Poisson-bracket techniques. They find the electron's trajectories condition- 
ally periodic^^ and regionally confined due to restrictions from E, M and flc- 
However, unlike Pauli, who used Sommerfeld quantization, SR employ the 
Einstein-Brillouin-Keller (EBK) quantization conditions. 



The background of EBK quantization touches on the foundations of clas- 
sical and quantum mechanics.'^ For the present purpose its essential rationale 
may be summarized as follows: A semiclassical treatment envolves turning 
points of radial, or other librating motion. Any tunneling through "forbid- 
den" regions of negative kinetic energy is ruled out. Viewed in terms of the 
quantum mechanical WKB (Wentzel-Kramers-Brillouin) approximation, the 
hard reflection of a wavefunction at a turning point corresponds to a so-called 
"loss" of phase (phase shift by vr) compared to the soft reflection caused by 
tunneling (phase shift by it/ 2 ). This shortcoming can be remedied with 
EBK quantization conditions by addition of a value of 1/4, for each libra- 
tional turning point, to the corresponding quantum number. Such is the case 
for the electron's elliptical and hyperbolic librations in the above quantiza- 
tion, Eq. (19ab), but not for a rotation about the z axis, Eq. (19c). Strand 




(19ab) 



and 




(19c) 
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and Reinhardt call these quantization conditions "primitive" to distinguish 
them from more sophisticated ones, specified below. 

A semiclassical treatment of a free H atom with EBK quantization has 
recently been presented in these pages.l^ In this case the isotropic symmetry 
permits a separation of variables in spherical coordinates, r, 9 and and 
the EBK quantization conditions are like Eqs. (19abc) except for j = r, 6. 
The atomic ground state is characterized by the quantum numbers (n^, ng, 
n^) = (0, 0, 0). Accordingly, the action in the atom's ground state, Ai = Ih, 
is attributed only to tunneling at the radial and latitudinal turning points 
(the former being the nucleus). Applying EBK quantization to the ground 
state of H2, denoted IsS^ in molecular spectroscopy, SR likewise assign the 
quantum numbers {n^,nri,n^) = (0,0,0). 

B. Results 

Energies of the Coulomb oscillator (CO), in adiabatic dependence on the 
proton separation, are listed in Appendix B and shown in Fig. 7 in com- 
parison with exact quantum mechanical (QM) results and the semiclassical 
calculation by SR.I^ The lower part of Fig. 7 shows the electronic energy 
Eei{R) of the ground state, IsEg. At large proton separations, R > 6 tb, 
both the SR calculation (circles) and the CO approach (crosses) agree ex- 
cellently with the exact QM values (curve). This is the situation where the 
electron stays near one nucleus (see Fig. 5). As Fig. 7 further shows, such 
agreement ceases once the classical electron motion leads beyond the molec- 
ular bisector, which happens for proton separations below the bifurcation 
value, R < R ^ 5.57 (see Fig. 6). Interestingly, the deviations of CO 
and SR from QM are opposite over the entire range. The SR results are 
discontinuous at a certain proton separation, R* ^ 1.38 r^. Remarkably, at 
(or near) that value, R*, the CO energy crosses the curve of the QM solution. 
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Fig. 7. Dependence of energies of on internuclear distance 
R. 

Bottom: Electronic energy £^ez from quantum mechanics (QM, 
solid curve), primitive semiclassical quantization by Strand and 
Reinhardt (SR, Q) ^'^'^ the present Coulomb-oscillator approach 
(CO, +). 

Middle: Total energy E of the IsS^ ground state from QM 
(dashed curve) with minimum (•), values by SR (□) and CO (x) 
and their average (o) , energy of a free H atom (dashed horizontal 
line). 

Top: Total energy E of the 2pY\.u state from QM (dotted curve) 
with minimum (•), and historical value by Pauli and Niessen (A). 



Adding to the electronic energy E^i the proton-proton repulsion gives the 
total energy Eq. (15). The middle part of Fig. 7 shows by the dashed 
curve the exact total energy E{R) of the IsE^ ground state, obtained from 
QM and, by symbols, the corresponding CO and SR values. The sohd dot at 
the minimum of the curve shows the QM equilibrium energy, in agreement 
with experiment, Eq — —1.20 and the equilibrium internuclear distance. 
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Rq = 2.00 Vb- The CO energy (x) comes out too low, due to the inaccuracy 
of its Eel, with an equihbrium value Eq{CO) ^ —1.38 Ry at Ro{CO) ^ 2.5 
tb- Conversely, the SR energy (□) comes out too high with Eq{SR) ~ —1.05 
Ry at Ro{SR) ~ 5 tb- Since the CO and SR results deviate about equal and 
oppositely from QM, their average {()) is close to the exact values with a 
minimum of So[|(CO + SR)] ^ -1.19 Ry at i?o[|(CO + SR)] ^ 2.5 r^. 

The molecular binding energy is the difference of E in the molecule and in 
the constituting atoms, here, Eq{H2)— E{H). The energy of a free hydrogen 
atom, E{H) = —1 Ry, is indicated in Fig. 7 by the fine horizontal line. 
Both semiclassical treatments, CO and (barely) SR, yield molecular binding 
energies with negative values and thus a stable molecule ion. Why are 
they more successful than the early attempts, in the 1920s, by Pauli,'^ and 
independently Niessen,'22l with the Sommerfeld quantization conditions of the 
old quantum theory? 

For the same reason that Sommerfeld had excluded the angular quan- 
tum number / = for the H atom — avoidance of electron collision with 
the nucleus — both Pauli and Niessen excluded electron motion in the nu- 
clear plane of the H2 molecule ion. They then found the lowest admissible 
quantum state to be (n^, n^, n^) = (0, 1, 1), denoted 2pn„ in molecular spec- 
troscopy, with Eo{P,N) = -0.52 Ry at Ro{P, N) = 5.53 ±0.01 r^, depicted 
by the triangle in the top part of Fig. 7. The QM energj^^U of that quantum 
state is Eo{2pIlu) = —0.27 Ry at Rq{2pIIu) ~ 8 r^, marked by the solid dot at 
the minimum of the dotted curve. Since both these energy values are higher 
than the ground state of a free hydrogen atom, they give rise to positive 
molecular binding energies and thus to spontaneous dissociation, H2{2pllu) 
^ H + . Qualitatively, Pauli's and Niessen's finding of energetic insta- 
bility is borne out by quantum mechanics for this exited state of H2 (Pauli's 
argumentP about "dynamical stability" notwithstanding). The deviation of 
their historical value (A) from the (dotted) QM curve is remarkably small — 
comparable to those of the CO and SR results for the ground state. Pauli's 
and Niessen's misfortune, though, was that they misinterpreted their result 
as the molecule ion's ground state — an assessment with fateful consequences 
in the development of quantum theory. 
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C. Discussion 



Why do the semiclassical results of the Coulomb oscillator and of SR's 
quantization deviate from the QM solution of the H2 molecule ion? Strand 
and Reinhardt explain the deviation of their primitive quantization from QM 
with effective potential barriers arising from constrictions due to conserva- 
tion of both the energy E and the bifocal Runge-Lenz component VLc (the 
angular momentum vanishes for the ground state, M = 0). The most drastic 
consequence of those barriers is the discontinuity of Eei at R* (see Fig. 7) 
and the large deviations at closer proton separation, R < R*. While the sim- 
ulation of quantum mechanical tunneling beyond the semiclassical turning 
points of librations is adequately achieved by EKB quantization under the 
isotropic symmetry of a free H atom^Fjq. (19) is less successful under the 
lower symmetry of H2 ■ When SR remedy the situation with "unified" semi- 
classical quantization conditions, then E^i agrees, for all practical purposes, 
with QM. Those unified quantization condition, going well beyond Eq. (19), 
are sophisticated in their dependence on E^i, flc, and the hyperbolic turning 
points 1]^. They will not be discussed here. 

The reason for the deviation of the CO results from the QM values is the 
neclect of the electron's wave nature in the underlying quantization condition, 
Eq. (16). In proposing his wave hypothesis de Brogli^^ already showed 
that the quantization condition of the Bohr model. An = nh, is equivalent 
to n standing waves along the nth Bohr orbit. If s denotes the position 
along the Bohr orbit, then the de Broglie wave can be expressed as w{s) = 
sm[2Tran{s)/h], with the variable an(s) = {An/ Sn) Jq ds' along the orbit's 
circumference Sn = '^nrn- A generalization gives the de Broglie wave of the 
Coulomb oscillator of a free H atom in the ground state (n = 1), 



where v{z') is the speed from Eq. (2) and the fold-out factor from Eqs. 



w{z) 



cos[27ia{z) / h] 



(20a) 



with 




(20b) 



(4), (9) and (14). 
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Fig. 8. de Broglie wave of the free H atom Coulomb oscillator 
(curve) and of the H2 Coulomb oscillator (small circles) for the 
same proton separation as in Fig. 5, R = 6 tb- Large circles 
indicate the axial positions of the nuclei. 

This de Broglie wave, shown by the line curve in Fig. 8 for a free H 
atom positioned at +c, has peak at the nucleus — like the QM radial wave 
function — and a node at each turning point. By the above characterization, 
those turning points arc "soft," and it is their softness that ensures the exact 
energy of the free atom. However, when Eq. (20ab) is applied to H2 for 
proton separations beyond the bifurcation value, R > R , with axial speed 
from Eq. (8ab) and integration away from the occupied nucleus, then 
the de Broglie wave is found to be "truncated" (no nodes) at both the outer 
and inner turning point (sec Fig. 8, small circles). Those turning points are 
"hard" and give rise to incorrect energies. Qualitatively, an augmentation 
of the truncated de Broglie wave with (exponential) "tunneling tails," deter- 
mined by the negative kinetic energy in the classically "forbidden" region. 
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would "soften" the turning points. This would give rise to an effective far 
turning point farther out, Cg// > C, and accordingly raise the CO energy, 
Eq. (7), toward the QM result. 




Fig. 9. de Broglie wave of the H2 Coulomb oscillator along axial 
distance OC from Fig. 4 (right side of graph) and lateral distance 
OB (left side). Large circles indicate the axial positions of the 
nuclei, here with the same separation as in Fig. 6, R — 2 rs- 

When the proton separation is below bifurcation, R < R, then Eq. (20b) 
should be integrated from the midpoint in both the axial and lateral direc- 
tion rather than favoring one nucleus with the crest of the de Broglie wave. 
Again, the de Broglie wave is found to be truncated at the axial and lateral 
turning points, C and B, respectively (see Fig. 9). An exception exits for 
the proton separation R*. The de Broglie wave, shown in Fig. 10, then has 
a minimum at both turning points, w{C) — w{B) — —1, according to action 
values of = and Ay = Such turning points seem to be "benign" — 
reminiscent of the soft turning points in the free-atom case — and cause the 
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CO energy Eei{R*) in Fig. 7 to be exact. At still smaller proton separation, 
R < R*. the de Broglie wave is truncated again (not shown). The limit R = 
corresponds to a Coulomb oscillator in the free He^ ion which, like in the 
free H atom for R = oo, has de Broglie nodes at the turning points and an 
exact energy value. 



w 




Fig. 10. Same as Fig. 9 but for the proton separation R* = 1.38 
rs where the semiclasssical energy is exact. 

If the explanation that the CO energy Egi{R) deviates from the QM 
curve because of the neglect of wave effects is valid, then this sheds new light 
on the SR results. The opposite sign of the CO and SR deviations then 
suggests that SR's primitive EBK quantization, while appropriate for a free 
atom, simulates too much wave effects under the lower symmetry of the H2 
molecule ion. However, the average of both those semiclassical quantizations 
seems to be an excellent compromise, as evidenced by the close agreement of 
the corresponding total energy (0) with the (dashed) QM curve in Fig. 7. 
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In conclusion, semiclassical quantization can rise to Gutzwiller's chal- 
lenge and "explain" the chemical bond in the paradigm molecule,'^ if^^, by a 
combination of classical mechanics, quantization, and moderate wave effects. 
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V. APPENDIX A: QUANTIZATION 

By Eq. (4) the action integral of the atomic Coulomb oscillator is 



A = \ (fpjz = -^V2Zm [ J- - ^dz. (4') 



r 



For a comparison with integral tables we change notation to x = z and use 
the abbreviation a = — l/P. Then 



11, /■ /I , f VX , 

az = / \ — \- aax = / ax (211 

^ r . \ X . X 



with X = ax"^ + X. Integration tables give 



The first term on the rhs, evaluated at the limits of the Jp integration, 
vanishes. The last integral in Eq. (22), tabulated as 

J = i-Vf) arcsin(l - 2a;/r), (23) 
and evaluated at the limits, x = and x = T, contributes 

— [arcsin(l) — arcsin(— 1)] a/T = —nVT. (24) 
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Combining Eqs. (4'), (22) and (24), together with a fold-out factor (j) — 2, 
gives the action integral, to be equated with the Sommerfeld quantization 
condition. 

An = -^V2Zm-{-7rVf) = nh. (25) 
(p 2 

We square Eq. (25) and solve for the amplitude of the quantized Coulomb 
oscillator, 

Tn = 2^n2, (26) 

in terms of the Bohr radius r^. 

VI. APPENDIX B: DATA 

TABLE I. Electronic ground-state energy Eei for various nuclear separations 
R of the H2 molecule from quantum-mechanical calculations (QM, Ref 21), 
the "primitive" semiclassical quantization of Strand and Reinhardt (SR, Ref. 
15), and the present Coulomb-oscillator approach (CO). 
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K [a.u.) 


QM (Ry) 


0.00 


-4.00 


0.25 


-3.80 


0.50 


-3.47 


0.75 


-3.11 


1.00 


-2.90 


1.25 


-2.68 


1.50 


-2.50 


1.75 


-2.34 


2.00 


-2.21 


2.50 


-1.99 


3.00 


-1.82 


3.50 


-1.69 


4.00 


-1.59 


5.00 


-1.45 


6.00 


-1.36 


8.00 


-1.26 


10.00 


-1.20 


12.00 


-1.17 



SR (Ry) 


CO {Ry 


-4.00 


-4.00 


-4.75 


-3.35 


-4.27 


-3.12 


-3.76 


-2.93 


-3.34 


-2.81 


-3.07 


-2.67 


-2.19 


-2.55 


-2.09 


-2.45 


-1.99 


-2.36 


-1.79 


-2.18 


-1.66 


-2.02 


-1.56 


-1.87 


-1.49 


-1.75 


-1.40 


-1.54 


-1.33 


-1.38 


-1.25 


-1.27 


-1.20 


-1.20 


-1.17 


-1.17 



23 



References 

^ Other failures were its inability to give the brightness of spectral lines and 
the unsuccessful extension to the He atom. 

^M. Jammer, The Conceptual Development of Quantum Mechanics 
(McGraw-Hill, New York, 1966), p. 129. 

^Instead of the traditional notation k for the angular quantum number in 
the old quantum theory, we will use the letter / (= k) to facilitate the 
connection with quantum mechanics. 

^W. Pauli, "Uber das Modell des Wasserstoffmolekiilions , " Annalen der 
Physik 68, 177-240 (1922). A history of Pauli's thesis is given by C. 
P. Enz, "No Time to he Brief: A Scientific Biography of Wolfgang Pauli" 
(Oxford UP, 2002), pp. 63-74. 

^M. C. Gutzwiller, Chaos in Classical and Quantum Mechanics (Springer, 
New York, 1990), p. 36. 

^A. Sommerfeld, "Zur Quantentheorie der Spektrallinien," Annalen der 
Physik, 51, 1-94 (1916). 

^This length of semiminor axis has been modified from Sommerfeld's original 
expression, bni = {rB/Z)nk, to achieve agreement with expressions from 
quantum mechanics. 

^If the nuclear charge is uniformly distributed, then Z'{r) — {r/r^YZ, and 
the electron's motion inside the nucleus is harmonic. 

^M. Bucher, D. Elm and D. P. Siemens, "Average position in Kepler mo- 
tion," Am. J. Phys. 66, 929-930 (1998). 

^°L. Pauling and E. B. Wilson, Introduction to Quantum Mechanics (Dover, 
New York, 1935) p. 144. 

^^Although almost all quantum texts use the weaker statement "at the nu- 
cleus." 

^^M. Bucher, "The electron inside the nucleus: An almost classical derivation 
of isotropic hyperfine interaction," Eur. J. Phys. 21, 19-22 (2000). 



24 



^^A. Einstein, "Zum Quantenansatz von Sommerfeld und Epstein," Verh. 
Dtsch. Phys. Ges. 19, 82-92 (1917). A summary with modem comments 
is given by A. D. Stone, "Einstein's unknown insight and the problem of 
quantizing chaos, " Phys. Today, Aug. 2005, pp. 37-43. 

^"^Eq. (3), introduced as Sommerfeld quantization for famiharity's sake, is 
also an Einstein quantization condition. 

^^M. P. Strand and W. P. Reinhardt, "Semiclassical quantization of the low 
lying electronic states of H^," J. Chem. Phys. 70, 3812-27 (1979). 

^^The (one-center) Runge-Lenz vector, f2 = v x L — Ze'^r/r, of a Sommer- 
feld orbit of angular momentum L has a magnitude proportional to the 
ellipse's eccentricity e and points from the nucleus toward the perihelion; 
see J. Morehead, ''Visualizing the extra symmetry of the Kepler prob- 
lem," Am. J. Phys. 73, 234-239 (2005). For a Coulomb oscillator (L = 0, 
e — 1), Vt oscillates and thus is not conserved. Neither is Qc of the two- 
center Coulomb oscillator conserved; see H. A. Erikson and E. L. Hill, 
"y4 Note on the one-electron state of diatomic molecules," Phys. Rev. 75, 
29-31 (1949); and also C. A. Coulson and A. Joseph, "^4 constant of the 
motion for the two-center Kepler problem, " Int. J. Quantum Chem. 1 , 
337-347 (1967). 

^^The trajectory will eventually revisit an arbitrarily close neighborhood of 
any previous point. 

^^Ref. 5, pp. 207-215. 

^^L. J. Curtis and D. G. Ellis, ''Use of the Einstein-Brillouin-Keller action 
quantization," Am. J. Phys. 72, 1521-1523 (2004). 

^°K. F. Niessen, "Zur Quantentheorie des Wasserstoffmolekulions, " Annalen 
der Physik 70, 129-134 (1923). 

2iM. M. Madsen and J. M. Peek, At. Data 2, 171 (1971); E. Teller and H. 
L. Sehlin, in Physical Chemistry, An Advanced Treatise (Academic, New 
York, 1970), Vol. 5, p. 35. 

22L. de Broghe, "Ondes et quanta," Comptes Rendus 177, 507-510 (1923). 



25 



K. Ruedenberg, "T/ie physical nature of the chemical bond," Rev. Mod. 
Phys. 34, 326-276 (1962). 



26 



2-24-2006 



Coulomb oscillation in the hydrogen atom 
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Abstract 

Semiclassical oscillation of the electron through the nucleus of the H atom 
yields both the exact energy and the correct orbital angular momentum for 
/ = quantum states. Similarly, electron oscillation through the nuclei of 
H2 accounts for a stable molecule ion with energy close to the quantum 
mechanical solution. The small discrepancy arises from the neglect of the 
electron's wave nature. 

PACS numbers: 03.65.Sq, 31.10.+Z, 31.20.Pv 
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I. INTRODUCTION 



Two of the reasons why the old quantum theory of Bohr and Sommerfeld 
was abandoned in the mid 1920s were the theory's failure to give the correct 
multiplet structure of the hydrogen atom and the stability of the hydrogen 
molecule ion, The old vector model of angular momentumP gave, for 

a given principal quantum number n, sublevels with angular quantum num- 
bers I = 1, 2, nP Spectroscopic evidence, however, showed multiplicities 
of spectral line splitting in a magnetic field according to / = 0, 1, ...,n — 1. 
Max Jammer, in his review,'^ notes that "the old quantum theory could never 
resolve this inconsistency." 

A treatment of the hydrogen molecule ion with Sommerfeld's quantization 
conditions had been Wolfgang Pauli's doctoral thesis of 1922.'^Pauli found its 
molecular binding energy to be positive (non-binding) — contrary to (later) 
experimental findings. Martin Gutzwillei'^ thinks that "the solution of this 
problem can be rated, with only slight exaggeration, as the most important 
in quantum mechanics, because if an energy level with a [more] negative 
value [than of a free hydrogen atom] can be found, then the chemical bond 
between two protons by a single electron has been explained." 

Both dilemmas of the old quantum theory can be resolved, though, with 
a single extension: an oscillation of the electron through the nucleus (nuclei) 
of the atom (molecule) . In essence this solution is already formally included 
in Sommerfeld's theory of the hydrogen atonP but was explicitly omitted by 
Sommerfeld and his school as being unphysical. The case in point, obtained 
with Sommerfeld's quantization conditions for radial and angular motion, 
is a quantum state with zero angular action, characterized by an angular 
quantum number / = 0. What is its orbit? 

The geometry of an nl Sommerfeld ellipse is given by its semimajor axis, 
dni = {j'B/Z)n^, and semiminor axis, hni = {rB/Z)n^^/l{l + 1)^ Here tb = 
/i^/47r^me^ is the Bohr radius in terms of fundamental constants, serving as 
an atomic distance unit, and Ze is the nuclear charge. An (n, 0) orbit is thus a 
line ellipse with its nuclear focus at one end and its empty focus at the other. 
This case was regarded as unphysical because of the electron's collision with 
the nucleus — an uncritical adaptation from celestrial mechanics. A closer 
inspection confirms that a line ellipse with terminal nuclear focus is indeed 
unphysical — but for quite a different reason! 

Figure 1 Here 
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Fig. 1. Partial trajectory of an extranuclear orbit XX' and of a penetrating 
orbit PP' through nucleus A^. The dotted line S shows the major symmetry axis 
of XX'. 

Leaving quantization conditions momentarily aside, what would happen 
if we continuously decrease an angular quantum number A while keeping the 
principal quantum number n constant? We then would get more and more 
slender ellipses with the same length of major axis, 2a„. By basic electric 
theory, the nuclear Coulomb potential outside a finite-size nucleus N is 
given by the point potential as if all nuclear charge, +Ze, was concentrated 
at the center of the nucleus. This holds as long as the electron orbit stays 
outside the nucleus. Two borderline cases are illustrated in Fig. 1. For a very 
small A value, say < Ax <^ 1, we obtain a very slender elliptical orbit with 
partial trajectory XX' about the nucleus. Further decrease of A to Ap < Ax 
causes an intrusion of the electron into the finite nucleus of radius tat (see 
trajectory PN in Fig. 1). Once inside, at a distance r < from the center, 
then, by Gauss's law, only a fraction of the nuclear charge, Z'{r)e < Ze, 
acts on the electron via centripetal force.l^ Accordingly, the electron's exit 
trajectory NP' is no longer symmetric to its approach trajectory PN with 
respect to the major axis S of the (partial) ellipse XX'. In the extreme 
case of a head-on penetration of the nucleus, A = 0, there is no centripetal 
force at all! The electron will then, with almost constant speed, traverse the 
nucleus, continue, with decreasing speed, to the opposite turning point of its 
line orbit and revert its motion periodically. We want to call the electron's 
straight-line oscillation in the Coulomb potential of a finite-size nucleus a 
"Coulomb oscillator. " 



II. HYDROGEN ATOM 

For the formal treatment of the non-relativistic Coulomb oscillator we 
designate the z axis along the line orbit, with turning points at z = ±r and 
nuclear position at z = (see Fig. 2). The electron's total energy E at 
position z must equal the potential energy at a turning point, 

E = -mv^ — = . (1) 

2 |2| r ^ ' 

This gives the electron's speed along the z axis. 
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^ = ±e\— Wo-T^- (2) 



Figure 3 displays its dependence on the axial position as a two-wing curve 
cusped at the nucleus. Atomic units (a.u.) are used, that is, the Bohr radius 
rs and the "Bohr speed" vb = 2'Ke^/h = ac — the electron speed in the 
ground-state Bohr orbit of the H atom — with fine-structure constant a ~ 
1/137 and speed of light c. The curve's wing along the positive z axis, 
l^l = r, gives the radial speed, Vr{r) — \v{z)\, necessary for Sommerfeld's 
radial quantization condition, 

Prdr — m <j) Vr{r)dr = rirh. (3) 

Here Pr is the radial momentum, rir — 1, 2, ... is the radial quantum number, 
and h is Planck's quantum of action. Integration is over one period of the 
radial motion, z = H-F — — > -|-F. Graphically, the radial quantization is 
illustrated in Fig. 3 by the area under one wing of the speed curve. For a 
line orbit, / = , the radial quantum number equals the principal quantum 
number, n = rir + I — rir- 

In order to express the radial quantization in terms of axial motion we 
employ a "fold-out factor," = 2, to compensate for the doubling of inte- 
gration range in the extension from the radial one-wing speed curve to the 
axial two-wing curve. The quantization is thus restated, 

1 f in f 

— f Pzdz = -J- f v{z)dz = Uzh, (3') 



0/" <f> 

with axial quantum number Uz = rir = n and integration over the axial 
double-wing range, z = +r — ^ — F — * +F. By symmetry we can restrict 
the axial action integral to one quarter of the oscillation, say z = +T — > 0, 



1 / , 4m /"V , , 4e /— — /l T, 

-- (t Pzdz — / v[z)dz — — -V2Zm / \ dz — nn. (4j 

9 J (p Jr <P 7rV^F 

Here the electron's motion in the negative z direction is accounted for by 
the negative sign. Although the electron's speed through a point nucleus 
diverges, v{0) — oo, the action integral, Eq. (4), stays finite and determines 
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the quantized amplitude r„ of the Coulomb oscillator. Graphically the am- 
plitude r„ must be such that it stretches the speed curve horizontally to 
the extent that the area under one wing, An = nh, represents the quantized 
action. The analytic solution, derived in Appendix A, is 

r. = 2^n2. (5) 
Inserting Eq. (5) into Eq. (1) yields the quantized energy. 



En = ^ = --rRy, (6) 

1 r) It' 



in terms of the Rydberg energy unit, Ry = 27T'^me^/h'^ = 13.6 eV, and in 
agreement with the energy of the nth Bohr orbit. 

Note that Eq. (5) gives the amplitude of the nth Coulomb oscillator 
as twice the radius of the nth Bohr orbit or of the semimajor axis of an 
nl Sommerfeld ellipse, r„ = = {rB/Z)n^. For comparison, the time- 
average radial distance of a Kepler orbiiP of major and minor semiaxes a 
and 6, respectively, is {r)t = (3a^ — 6^)/(2a). A line ellipse (6 = 0) has then 
{f)t = §«• Applied to an nl Sommerfeld orbit,'^ its average size, (r„;)j = 
(rB/Z)[3n^ — /(/ + l)]/2, is in agreement with the corresponding quantity 
from quantum mechanics,'^ (r„;) = / %l)*ripd?r. Thus the time-average radial 
distance of a Coulomb oscillator is (r^o)* = b / Z)n'^ . For the ground state 
of the hydrogen atom, n = 1, this gives (rio)t = Ir^, as is well-known from 
quantum mechanics.'^ 

The concept of the electron's semiclassical Coulomb oscillation is consis- 
tent with the Fermi-contact term of hyperfine interaction for / = states, 
which arises from the presence of the electron inside the nucleus. This is 
familiar from quantum mechanic^ and can be interpreted semiclassically as 
a local- field effect.'^ 

To be sure, the extension of Sommerfeld's theory by the Coulomb oscil- 
lator resolves the discrepancy of the old quantum theory with spectroscopy, 
mentioned above, only at the low end of angular quantum numbers, / = 0. 
The resolution at the high end — repeal of the circular Bohr orbit, / 7^ n — 
involves an analysis in terms of space quantization which is beyond the scope 
of the present study. 
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III. HYDROGEN MOLECULE ION 



A. Formalism 

A hydrogen molecule ion, H2, consists of two proton nuclei and one 
electron. We assume, in adiabatic approximation, the protons located at fixed 
positions z = ±c on the molecular axis (see Fig. 4). The term "Coulomb 
oscillator" denominates again the motion of a point electron, now along either 
the line through the protons {z axis) and with turning points ±C, or along 
the perpendicular line through the midpoint (y axis) with turning points ±B. 
The molecule ion's total energy E at any position on the axis, — C ^ ^ +C, 
must equal the potential energy at the turning point C, 



E = -mv^ - - r H = — \ . (7) 

2 \z + c\ \z-c\ 2c C + c C-c 2c ^ ^ 

For electron positions beyond the protons, z > c, this gives an electron speed 



, 2e / ^ C 



For positions between the protons, < 2; < c, the corresponding speed is 



, 2e / c C 



The speed expressions will be used in the action integral, 
1 /■ , 2m 





with outer contribution 



2m 

■^out — ~ 



and inner contribution 



/2m f^'^ 
p^dz = y v{z)dz = = Aout + Ain, (9) 



j Vout{z)dz + J Vout{z)dz 



(10a) 



An = I Vin{z)dz. (10b) 
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Here is a fold-out factor to be specified below. 

If the conditions are such that the electron swings along the z axis through 
the midpoint 0, that is, f j„(0) > 0, then there exists also an oscillation along 
the y axis with lateral speed m, having the same total energy, 

E = -mu^ H = H . (11) 

2 v^2^ 2c y/WT^ 2c ^ ^ 

Equal energy at the axial and lateral turning points, E{C) = E{B), Eqs. (7) 
and (11), determines the latters' geometric dependence. 



c^ 



5= + — -3c2. (12) 



Solving Eq. (11) for the lateral speed. 



_ 2e / 1 1 

^7^\l V7T^ ~ 7WT^' ^^^^ 

provides the integrand of the action integral over a lateral oscillation, 



2m 



-B 



B 



Subtraction of the protons' mutual repulsion from the total energy E of 
axial or lateral motion, Eqs. (7) and (11), gives the electronic energy. 



e2 



Eei = E- -, (15) 

in its dependence on the proton separation, R = 2c. This brackets the 
molecular problem with known atomic results, Eq. (6), in the limits of i? = cxo 
(free H atom, Z = 1) and R = (free He~^ ion, Z = 2). For those cases, as 
well as any proton-proton distance R between, we keep the action constant, 

A = + Ay = nh. (16) 

Equation (16) is the Einstein quantization conditiorP^' — a generalization of 
Sommerfeld's quantization over separable variables — where the quantum sum 
equals the sum of action integrals over topologically independent paths in 
phase space 
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Here we treat the molecule ion only in its ground state, n — 1. Analytic 
solutions of the action integrals. Eqs. (lOab) and (14), are complicated due 
to elliptic functions. We therefore integrate numerically and visualize the 
integrals by the area under the corresponding speed curves. The bold curve 
in Fig. 5 shows the axial electron speed v{z) for a far proton separation, 
R = 6rB- The electron, in its semiclassical motion, then oscillates only about 
(and through) the right proton. The area under the speed curve, Eq. (Sab), 
proportionally represents the ground state's unity of action, m ^ v{z)dz = (f)h, 
with a fold-out factor = 2 in analogy to the free-atom case, Eq. (4). The 
thin curve shows, for comparison, the axial electron speed in a free H atom — 
familiar from Fig. 3 — centered at the same proton position, +c. The pull 
from the left proton (at — c) on the oscillating electron can be seen by the 
distortion of the speed function v{z) and the redistribution of the area under 
the curve. 

When, with closer proton separation R, as in Fig. 6, the electron swings 
past the midpoint 0, then the single-cusp speed curve v{z) from Fig. 5 — 
akin in shape to letter A — becomes double-cusped (akin to letter M) and 
symmetric with respect to the bisector {y axis). Now there is also a lateral 
oscillation with speed u{y) having the same total energy E. The equality of 
E in both cases can be seen in Fig. 6 by the equahty of axial and lateral speed 
at the midpoint, ^;(0) = u{0) — a position where the electron experiences the 
same potential in either case. For convenience the lateral speed u{y), though 
perpendicular to the proton axis, is displayed in Fig. 6 together with the axial 
speed v{z). The lateral speed curve u{y), drawn dotted, is readily recognized 
by its dome shape (fl). At the bifurcation value of the proton separation, R, 
where the electron starts swinging though the midpoint 0, the axial speed 
curve v{z) changes from its one-centered A shape to a two-centered M shape. 
The area under the axial speed curve then abruptly doubles, M(i? — 6) ^ 
2A{R + d), upon a very small change in proton separation, 6 R. In order 
to keep the action integral continuous at R, the sudden area doubling is 
compensated by a corresponding doubling of the fold-out factor from (f) — 2 
for R > R to (j) = 4 for R < R. The same fold-out factor, = 4, must be 
used for the lateral action integral Ay, Eq. (14), as will become clear shortly. 

The unity of action, A = Ih, is visualized again in Fig. 6. To this end we 
compare the right half of the bold M-shape curve of axial electron speed in 
the molecule with the thin curve A (if) of the axial speed in a free if atom 
positioned at the right nucleus, +c. Due to attraction from the left nucleus. 
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the right wing of the M curve is smaller than that of A(iJ) by the area of 
lobe L. On the other hand, the left flank of the A{H) curve that extends over 
the negative z axis is smaller than the left wing of the lateral speed curve fl 
by the area of slice S. The area under the free-atom curve is then 

A{H) ^ Im + L + i n -5. (17a) 

The areas of lobe and slice are comparable, 

L^S. (17b) 

When the tiny notch to the right of the saddle point of M is taken into 
account, then the approximations (17ab) become equations and combine to 

M + n = 2A(if). (18) 

The area under both the axial and lateral speed curves is thus four times 
the area under one wing of the free-atom curve, M+fl = 4 x |A(if). Since 
the latter represents one quantum of action, h, the combined area M+fl 
visualizes its double fold-out, = 4. 

With very close proximity of the nuclei, R ^ 0, the crests of the M curve 
start merging while its saddle point, f (0), keeps rising. In the R = limit of 
fusing nuclei the axial electron speed becomes that of a free He~^ ion, M— »• 
A{He~^), familiar from Fig. 3. Concurrently, the lateral speed curve fl rises 
at its peak, u{0) = f (0), and narrows at its base until it, too, turns into the 
speed curve of the free He~^ ion, fl —>■ A{He~^). In the R = limit the three 
curves merge, M(0) = n(0) = A(iJe+). 

The results of the Coulomb-oscillator approach will be compared with 
another semiclassical calculation of by Strand and Reinhardt.!^ These 
authors, like Pauli,'^ separate the equation of motion in spheroidal coordi- 
nates, ^ = (r_|_ + r_)/2c , T] = (r+ — r_)/2c and ip by virtue of the constants 
of the motion: total energy E, angular momentum M about the z axis, and 
a component of the bifocal Runge-Lenz vector, ^Ic^ Here r+ (r_) is the 
distance of the electron from the nucleus at +c (— c). Strand and Reinhardt 
(SR) solve the ensuing one-dimensional differential equations with classical 
Poisson-bracket techniques. They find the electron's trajectories condition- 
ally periodic!^ and regionally confined due to restrictions from E, M and fic- 
However, unlike Pauli, who used Sommerfeld quantization, SR employ the 
Einstein-Brillouin-Keller (EBK) quantization conditions. 
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Pjdj = (uj + -)h, j = C, V 



(19ab) 



and 



n^h. 



(19c) 



The background of EBK quantization touches on the foundations of clas- 
sical and quantum mechanics.^ For the present purpose its essential rationale 
may be summarized as follows: A semiclassical treatment envolves turning 
points of radial, or other librating motion. Any tunneling through "forbid- 
den" regions of negative kinetic energy is ruled out. Viewed in terms of the 
quantum mechanical WKB (Wentzel-Kramers-Brillouin) approximation, the 
hard reflection of a wavefunction at a turning point corresponds to a so-called 
"loss" of phase (phase shift by vr) compared to the soft reflection caused by 
tunneling (phase shift by 7r/2 ). This shortcoming can be remedied with 
EBK quantization conditions by addition of a value of 1/4, for each libra- 
tional turning point, to the corresponding quantum number. Such is the case 
for the electron's elliptical and hyperbolic librations in the above quantiza- 
tion, Eq. (19ab), but not for a rotation about the z axis, Eq. (19c). Strand 
and Reinhardt call these quantization conditions "primitive" to distinguish 
them from more sophisticated ones, specified below. 

A semiclassical treatment of a free H atom with EBK quantization has 
recently been presented in these pages.l^ In this case the isotropic symmetry 
permits a separation of variables in spherical coordinates, r, 9 and y?, and 
the EBK quantization conditions are like Eqs. (19abc) except for j = r, 9. 
The atomic ground state is characterized by the quantum numbers (n^, n^, 
fT-if) = (0, 0, 0). Accordingly, the action in the atom's ground state, Ai = Ih, 
is attributed only to tunneling at the radial and latitudinal turning points 
(the former being the nucleus). Applying EBK quantization to the ground 
state of H2, denoted lsJ2g in molecular spectroscopy, SR likewise assign the 
quantum numbers {n^,nri,n^) = (0,0,0). 

B. Results 

Energies of the Coulomb oscillator (CO), in adiabatic dependence on the 
proton separation, are listed in Appendix B and shown in Fig. 7 in com- 
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parison with exact quantum mechanical (QM) results and the semiclassical 
calculation by SR.!^ The lower part of Fig. 7 shows the electronic energy 
Eei{R) of the ground state, IsEg. At large proton separations, R > 6 r^, 
both the SR calculation (circles) and the CO approach (crosses) agree ex- 
cellently with the exact QM values (curve). This is the situation where the 
electron stays near one nucleus (see Fig. 5). As Fig. 7 further shows, such 
agreement ceases once the classical electron motion leads beyond the molec- 
ular bisector, which happens for proton separations below the bifurcation 
value, R < R ^ 5.57 rs (see Fig. 6). Interestingly, the deviations of CO 
and SR from QM are opposite over the entire range. The SR results are 
discontinuous at a certain proton separation, R* ^ 1.38 r^. Remarkably, at 
(or near) that value, R*, the CO energy crosses the curve of the QM solution. 

Adding to the electronic energy Egi the proton-proton repulsion gives the 
total energy E, Eq. (15). The middle part of Fig. 7 shows by the dashed 
curve the exact total energy E{R) of the IsS^ ground state, obtained from 
QM and, by symbols, the corresponding CO and SR values. The solid dot at 
the minimum of the curve shows the QM equilibrium energy, in agreement 
with experiment, Eq = —1.20 Ry, and the equilibrium internuclear distance, 
Ro = 2.00 tb- The CO energy (x) comes out too low, due to the inaccuracy 
of its Eel, with an equilibrium value Eq{CO) ~ —1.38 Ry at Ro{CO) ~ 2.5 
tb- Conversely, the SR energy (□) comes out too high with Eq{SR) ~ —1.05 
Ry at Ro{SR) ~ 5 r^. Since the CO and SR results deviate about equal and 
oppositely from QM, their average {(}) is close to the exact values with a 
minimum of ^o[|(CO + SR)] « -1.19 Ry at i?o[|(CO + SR)] ^ 2.5 r^. 

The molecular binding energy is the difference of E in the molecule and in 
the constituting atoms, here, EqIH^)— E{H). The energy of a free hydrogen 
atom, E{H) = —1 Ry, is indicated in Fig. 7 by the fine horizontal line. 
Both semiclassical treatments, CO and (barely) SR, yield molecular binding 
energies with negative values and thus a stable molecule ion. Why are 
they more successful than the early attempts, in the 1920s, by Pauli,^ and 
independently Niessen,'^with the Sommerfeld quantization conditions of the 
old quantum theory? 

For the same reason that Sommerfeld had excluded the angular quan- 
tum number / = for the H atom — avoidance of electron collision with 
the nucleus — both Pauli and Niessen excluded electron motion in the nu- 
clear plane of the H2 molecule ion. They then found the lowest admissible 
quantum state to be (n^, n^, n^) = (0, 1, 1), denoted 2pllu in molecular spec- 
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troscopy, with Eo{P,N) = -0.52 Ry at Ro{P,N) = 5.53 ±0.01 r^, depicted 
by the triangle in the top part of Fig. 7. The QM energy^ of that quantum 
state is Eo{2pIlu) = —0.27 Ry at Rq{2pIIu) ~ 8 r^, marked by the sohd dot at 
the minimum of the dotted curve. Since both these energy values are higher 
than the ground state of a free hydrogen atom, they give rise to positive 
molecular binding energies and thus to spontaneous dissociation, if^(2pnu) 
^ H + H^. Qualitatively, Pauli's and Niessen's finding of energetic insta- 
bility is borne out by quantum mechanics for this exited state of H2 (Pauli's 
argument^ about "dynamical stability" notwithstanding). The deviation of 
their historical value (A) from the (dotted) QM curve is remarkably small — 
comparable to those of the CO and SR results for the ground state. Pauli's 
and Niessen's misfortune, though, was that they misinterpreted their result 
as the molecule ion's ground state — an assessment with fateful consequences 
in the development of quantum theory. 

C. Discussion 

Why do the semiclassical results of the Coulomb oscillator and of SR's 
quantization deviate from the QM solution of the H2 molecule ion? Strand 
and Reinhardt explain the deviation of their primitive quantization from QM 
with effective potential barriers arising from constrictions due to conserva- 
tion of both the energy E and the bifocal Runge-Lenz component flc (the 
angular momentum vanishes for the ground state, M = 0). The most drastic 
consequence of those barriers is the discontinuity of E^i at R* (see Fig. 7) 
and the large deviations at closer proton separation, R < R*. While the sim- 
ulation of quantum mechanical tunneling beyond the semiclassical turning 
points of librations is adequately achieved by EKB quantization under the 
isotropic symmetry of a free H atom^ Eq. (19) is less successful under the 
lower symmetry of H2 ■ When SR remedy the situation with "unified" semi- 
classical quantization conditions, then E^i agrees, for all practical purposes, 
with QM. Those unified quantization condition, going well beyond Eq. (19), 
are sophisticated in their dependence on Egi, Qc, and the hyperbolic turning 
points ?7_|_. They will not be discussed here. 

The reason for the deviation of the CO results from the QM values is the 
neclect of the electron's wave nature in the underlying quantization condition, 
Eq. (16). In proposing his wave hypothesis de Broglie^^ already showed 
that the quantization condition of the Bohr model. An = nh, is equivalent 
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to n standing waves along the nth Bohr orbit. If s denotes the position 
along the Bohr orbit, then the de Broglie wave can be expressed as w{s) = 
sm[27ian{s)/h], with the variable an(s) = (A„/S'„) Jq ds' along the orbit's 
circumference Sn = STrr^. A generalization gives for the Coulomb oscillator 
of a free H atom in the ground state {n = 1), its de Broglie wave as 



and the speed v{z') from Eq. (2). This de Broglie wave, shown by the line 
curve in Fig. 8 for a free H atom positioned at +c, has a node at the nucleus 
and at each turning point. By the above characterization, those turning 
points are "soft," and it is their softness that ensures the exact energy of the 
free atom. However, when Eq. (20ab) is applied to H2 for proton separations 
beyond the bifurcation value, R > R, with axial speed from Eq. (Sab) 
and integration away from the occupied nucleus, then the de Broglie 

wave is found to be "truncated" (no nodes) at both the outer and inner 
turning point (see Fig. 8, small circles). Those turning points are "hard" 
and give rise to incorrect energies. Qualitatively, an augmentation of the 
truncated de Broglie wave with (exponential) "tunneling tails," determined 
by the negative kinetic energy in the classically "forbidden" region, would 
"soften" the turning points. This would give rise to an effective far turning 
point farther out, Ce// > C, and accordingly raise the CO energy, Eq. (7), 
toward the QM result. 

When the proton separation is below bifurcation, R < R, then Eq. (20b) 
should be integrated from the midpoint in both the axial and lateral direc- 
tion rather than favoring one nucleus with the node of the de Broglie wave. 
Again, the de Broglie wave is found to be truncated at the axial and lateral 
turning points, C and B, respectively (see Fig. 9). An exception exits for 
the proton separation R*. The de Broglie wave, shown in Fig. 10, then has 
a minimum at both turning points, w{C) = w{B) = —1, according to action 
values of = |/i and Ay = ^h. Such turning points seem to be "benign" — 
reminiscent of the soft turning points in the free-atom case — and cause the 
CO energy Eei{R*) in Fig. 7 to be exact. At still smaller proton separation, 
R < R*, the de Broglie wave is truncated again (not shown). The limit R = 
corresponds to a Coulomb oscillator in the free He~^ ion which, like in the 



w{z) 



sin[27ra(^)//i] 



(20a) 



with 




(20b) 
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free H atom for R = oo, has de Broglie nodes at the turning points and an 
exact energy value. 

If the explanation that the CO energy Eei{R) deviates from the QM 
curve because of the neglect of wave effects is valid, then this sheds new light 
on the SR results. The opposite sign of the CO and SR deviations then 
suggests that SR's primitive EBK quantization, while appropriate for a free 
atom, simulates too much wave effects under the lower symmetry of the H2 
molecule ion. However, the average of both those semiclassical quantizations 
seems to be an excellent compromise, as evidenced by the close agreement of 
the corresponding total energy (0) with the (dashed) QM curve in Fig. 7. 

In conclusion, semiclassical quantization can rise to Gutzwiller's chal- 
lenge and "explain" the chemical bond in the paradigm molecule,'^ ifg^, by a 
combination of classical mechanics, quantization, and moderate wave effects. 



I thank Duane Siemens and Ernst Mohler for valuable discussions. Many 
thanks to Preston Jones for help with computer integration and graphics. I 
also thank Professor Gutzwiller for advice and his kind encouragement. 



By Eq. (4) the action integral of the atomic Coulomb oscillator is 




For a comparison with integral tables we change notation to x = z and use 
the abbreviation a = — l/F. Then 
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Vx 



X 



dx 



(21) 



with X = ax^ + X. Integration tables give 




(22) 
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The first term on the rhs, evaluated at the hmits of the Jp integration, 
vanishes. The last integral in Eq. (22), tabulated as 



/ 



'^'^ {-Vf) arcsin(l - 2x/r), (23) 



and evaluated at the limits, x — and x — F, contributes 

— [arcsin(l) — arcsin(— l)]v^=-7rv^. (24) 

Combining Eqs. (4'), (22) and (24), together with a fold-out factor (f) — 2, 
gives the action integral, to be equated with the Sommerfeld quantization 
condition, 

^„ = -^\/2Zm-(-7r\/f) = nh. (25) 
(p 2 

We square Eq. (25) and solve for the amplitude of the quantized Coulomb 
oscillator. 



r„ = 2^r^^ (26) 



in terms of the Bohr radius r^. 



VI. APPENDIX B: DATA 



TABLE I. Electronic ground-state energy E^i for various nuclear separations 
R of the H2 molecule from quantum-mechanical calculations (QM, Rcf 21), 
the "primitive" semiclassical quantization of Strand and Reinhardt (SR, Ref. 
15), and the present Coulomb-oscillator approach (CO). 
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K [a.u.) 


QM (Ry) 


0.00 


-4.00 


0.25 


-3.80 


0.50 


-3.47 


0.75 


-3.11 


1.00 


-2.90 


1.25 


-2.68 


1.50 


-2.50 


1.75 


-2.34 


2.00 


-2.21 


2.50 


-1.99 


3.00 


-1.82 


3.50 


-1.69 


4.00 


-1.59 


5.00 


-1.45 


6.00 


-1.36 


8.00 


-1.26 


10.00 


-1.20 


12.00 


-1.17 



SR (Ry) 


CO {Ry 


-4.00 


-4.00 


-4.75 


-3.35 


-4.27 


-3.12 


-3.76 


-2.93 


-3.34 


-2.81 


-3.07 


-2.67 


-2.19 


-2.55 


-2.09 


-2.45 


-1.99 


-2.36 


-1.79 


-2.18 


-1.66 


-2.02 


-1.56 


-1.87 


-1.49 


-1.75 


-1.40 


-1.54 


-1.33 


-1.38 


-1.25 


-1.27 


-1.20 


-1.20 


-1.17 


-1.17 
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VII. FIGURE CAPTIONS 



Fig. 1. Partial trajectory of an extranuclear orbit XX' and of a pene- 
trating orbit PP' through nucleus N . The dotted line S shows the major 
symmetry axis of XX' . 

Fig. 2. Line orbit of a Coulomb oscillator with nucleus at origin and 
turning points at ±r. 

Fig. 3. Axial electron speed v vs. position z of an electron in Coulomb 
oscillation for the ground state, n = 1, of a hydrogen atom H {Z = 1, solid 
curve) and a helium ion He'^ {Z — 2, dashed curve). The area under one 
wing of each curve represents the radial action, A—lh. 

Fig. 4. Axial Coulomb oscillation in an H2 molecule ion between ax- 
ial turning points ±C and through nuclei at ±c; perpendicular oscillation 
between lateral turning points ±B and through midpoint 0. 

Fig. 5. Axial speed v vs. position z in the ground state of an H2 molecule 
ion (bold curve) and, for comparison, of a free H atom (thin curve). Circles 
indicate the axial the positions of the nuclei, here with a large separation. 

Fig. 6. Axial speed v vs. position z in the ground state of an H2 
molecule ion (M-shaped bold curve) and, for comparison, of a free H atom 
(A-shaped thin curve, centered at the right nucleus). The Pi -shaped dotted 
curve, centered at 0, shows the lateral speed u vs. the perpendicular position 
y in the molecule. Circles indicate the axial positions of the nuclei, here with 
a small separation, R = 2 rs- 

Fig. 7. Dependence of energies of H2 on internuclear distance R. Bottom: 
Electronic energy E^i from quantum mechanics (QM, solid curve), primitive 
semiclassical quantization by Strand and Reinhardt (SR, Q) ^^'^ ^^e present 
Coulomb-oscillator approach (CO, -|-). Middle: Total energy E of the IsEg 
ground state from QM (dashed curve) with minimum (•), values by SR (□) 
and CO (x) and their average (o), energy of a free H atom (dashed horizontal 
line). Top: Total energy E of the 2pllu state from QM (dotted curve) with 
minimum (•), and historical value by Pauli and Niessen (A). 

Fig. 8. de Broghe wave of the free H atom Coulomb oscillator (curve) and 
of the Coulomb oscillator (small circles) for the same proton separation 
as in Fig. 5, R = 6 tb- Large circles indicate the axial positions of the nuclei. 

Fig. 9. de Broglie wave of the H2 Coulomb oscillator along axial distance 
OC from Fig. 4 (right side of graph) and lateral distance OB (left side). 
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Large circles indicate the axial positions of the nuclei, here with the same 
separation as in Fig. 6, i? = 2 r^. 

Fig. 10. Same as Fig. 9 but for the proton separation R* — 1.38 
where the semiclasssical energy is exact. 
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Abstract 

Semiclassical oscillation of the electron through the nucleus of the H atom 
yields both the exact energy and the correct orbital angular momentum for 
/ = quantum states. Similarly, electron oscillation through the nuclei of 
H2 accounts for a stable molecule ion with energy close to the quantum 
mechanical solution. The small discrepancy arises from the neglect of the 
electron's wave nature. 

PACS numbers: 03.65.Sq, 31.10.+Z, 31.20.Pv 
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I. INTRODUCTION 



Two of the reasons why the old quantum theory of Bohr and Sommerfeld 
was abandoned in the mid 1920s were the theory's failure to give the correct 
multiplet structure of the hydrogen atom and the stability of the hydrogen 
molecule ion, The old vector model of angular momentumP gave, for 

a given principal quantum number n, sublevels with angular quantum num- 
bers I = 1, 2, nP Spectroscopic evidence, however, showed multiplicities 
of spectral line splitting in a magnetic field according to / = 0, 1, ...,n — 1. 
Max Jammer, in his review,'^ notes that "the old quantum theory could never 
resolve this inconsistency." 

A treatment of the hydrogen molecule ion with Sommerfeld's quantization 
conditions had been Wolfgang Pauli's doctoral thesis of 1922.'^Pauli found its 
molecular binding energy to be positive (non-binding) — contrary to (later) 
experimental findings. Martin Gutzwillei'^ thinks that "the solution of this 
problem can be rated, with only slight exaggeration, as the most important 
in quantum mechanics, because if an energy level with a [more] negative 
value [than of a free hydrogen atom] can be found, then the chemical bond 
between two protons by a single electron has been explained." 

Both dilemmas of the old quantum theory can be resolved, though, with 
a single extension: an oscillation of the electron through the nucleus (nuclei) 
of the atom (molecule) . In essence this solution is already formally included 
in Sommerfeld's theory of the hydrogen atonP but was explicitly omitted by 
Sommerfeld and his school as being unphysical. The case in point, obtained 
with Sommerfeld's quantization conditions for radial and angular motion, 
is a quantum state with zero angular action, characterized by an angular 
quantum number / = 0. What is its orbit? 

The geometry of an nl Sommerfeld ellipse is given by its semimajor axis, 
dni = {j'B/Z)n^, and semiminor axis, hni = {rB/Z)n^^/l{l + 1)^ Here tb = 
/i^/47r^me^ is the Bohr radius in terms of fundamental constants, serving as 
an atomic distance unit, and Ze is the nuclear charge. An (n, 0) orbit is thus a 
line ellipse with its nuclear focus at one end and its empty focus at the other. 
This case was regarded as unphysical because of the electron's collision with 
the nucleus — an uncritical adaptation from celestrial mechanics. A closer 
inspection confirms that a line ellipse with terminal nuclear focus is indeed 
unphysical — but for quite a different reason! 

Figure 1 Here 



2 



Fig. 1. Partial trajectory of an extranuclear orbit XX' and of a penetrating 
orbit PP' through nucleus A^. The dotted line S shows the major symmetry axis 
of XX'. 

Leaving quantization conditions momentarily aside, what would happen 
if we continuously decrease an angular quantum number A while keeping the 
principal quantum number n constant? We then would get more and more 
slender ellipses with the same length of major axis, 2a„. By basic electric 
theory, the nuclear Coulomb potential outside a finite-size nucleus N is 
given by the point potential as if all nuclear charge, +Ze, was concentrated 
at the center of the nucleus. This holds as long as the electron orbit stays 
outside the nucleus. Two borderline cases are illustrated in Fig. 1. For a very 
small A value, say < Ax <^ 1, we obtain a very slender elliptical orbit with 
partial trajectory XX' about the nucleus. Further decrease of A to Ap < Ax 
causes an intrusion of the electron into the finite nucleus of radius tat (see 
trajectory PN in Fig. 1). Once inside, at a distance r < from the center, 
then, by Gauss's law, only a fraction of the nuclear charge, Z'{r)e < Ze, 
acts on the electron via centripetal force.l^ Accordingly, the electron's exit 
trajectory NP' is no longer symmetric to its approach trajectory PN with 
respect to the major axis S of the (partial) ellipse XX'. In the extreme 
case of a head-on penetration of the nucleus, A = 0, there is no centripetal 
force at all! The electron will then, with almost constant speed, traverse the 
nucleus, continue, with decreasing speed, to the opposite turning point of its 
line orbit and revert its motion periodically. We want to call the electron's 
straight-line oscillation in the Coulomb potential of a finite-size nucleus a 
"Coulomb oscillator. " 



II. HYDROGEN ATOM 

For the formal treatment of the non-relativistic Coulomb oscillator we 
designate the z axis along the line orbit, with turning points at z = ±r and 
nuclear position at z = (see Fig. 2). The electron's total energy E at 
position z must equal the potential energy at a turning point, 

E = -mv^ — = . (1) 

2 |2| r ^ ' 

This gives the electron's speed along the z axis. 
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^ = ±e\— Wo-T^- (2) 



Figure 3 displays its dependence on the axial position as a two-wing curve 
cusped at the nucleus. Atomic units (a.u.) are used, that is, the Bohr radius 
rs and the "Bohr speed" vb = 2'Ke^/h = ac — the electron speed in the 
ground-state Bohr orbit of the H atom — with fine-structure constant a ~ 
1/137 and speed of light c. The curve's wing along the positive z axis, 
l^l = r, gives the radial speed, Vr{r) — \v{z)\, necessary for Sommerfeld's 
radial quantization condition, 

Prdr — m <j) Vr{r)dr = rirh. (3) 

Here Pr is the radial momentum, rir — 1, 2, ... is the radial quantum number, 
and h is Planck's quantum of action. Integration is over one period of the 
radial motion, z = H-F — — > -|-F. Graphically, the radial quantization is 
illustrated in Fig. 3 by the area under one wing of the speed curve. For a 
line orbit, / = , the radial quantum number equals the principal quantum 
number, n = rir + I — rir- 

In order to express the radial quantization in terms of axial motion we 
employ a "fold-out factor," = 2, to compensate for the doubling of inte- 
gration range in the extension from the radial one-wing speed curve to the 
axial two-wing curve. The quantization is thus restated, 

1 f in f 

— f Pzdz = -J- f v{z)dz = Uzh, (3') 



0/" <f> 

with axial quantum number Uz = rir = n and integration over the axial 
double-wing range, z = +r — ^ — F — * +F. By symmetry we can restrict 
the axial action integral to one quarter of the oscillation, say z = +T — > 0, 



1 / , 4m /"V , , 4e /— — /l T, 

-- (t Pzdz — / v[z)dz — — -V2Zm / \ dz — nn. (4j 

9 J (p Jr <P 7rV^F 

Here the electron's motion in the negative z direction is accounted for by 
the negative sign. Although the electron's speed through a point nucleus 
diverges, v{0) — oo, the action integral, Eq. (4), stays finite and determines 
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the quantized amplitude r„ of the Coulomb oscillator. Graphically the am- 
plitude r„ must be such that it stretches the speed curve horizontally to 
the extent that the area under one wing, An = nh, represents the quantized 
action. The analytic solution, derived in Appendix A, is 

r. = 2^n2. (5) 
Inserting Eq. (5) into Eq. (1) yields the quantized energy. 



En = ^ = --rRy, (6) 

1 r) It' 



in terms of the Rydberg energy unit, Ry = 27T'^me^/h'^ = 13.6 eV, and in 
agreement with the energy of the nth Bohr orbit. 

Note that Eq. (5) gives the amplitude of the nth Coulomb oscillator 
as twice the radius of the nth Bohr orbit or of the semimajor axis of an 
nl Sommerfeld ellipse, r„ = = {rB/Z)n^. For comparison, the time- 
average radial distance of a Kepler orbiiP of major and minor semiaxes a 
and 6, respectively, is {r)t = (3a^ — 6^)/(2a). A line ellipse (6 = 0) has then 
{f)t = §«• Applied to an nl Sommerfeld orbit,'^ its average size, (r„;)j = 
(rB/Z)[3n^ — /(/ + l)]/2, is in agreement with the corresponding quantity 
from quantum mechanics,'^ (r„;) = / %l)*ripd?r. Thus the time-average radial 
distance of a Coulomb oscillator is (r^o)* = b / Z)n'^ . For the ground state 
of the hydrogen atom, n = 1, this gives (rio)t = Ir^, as is well-known from 
quantum mechanics.'^ 

The concept of the electron's semiclassical Coulomb oscillation is consis- 
tent with the Fermi-contact term of hyperfine interaction for / = states, 
which arises from the presence of the electron inside the nucleus. This is 
familiar from quantum mechanic^ and can be interpreted semiclassically as 
a local- field effect.'^ 

To be sure, the extension of Sommerfeld's theory by the Coulomb oscil- 
lator resolves the discrepancy of the old quantum theory with spectroscopy, 
mentioned above, only at the low end of angular quantum numbers, / = 0. 
The resolution at the high end — repeal of the circular Bohr orbit, / 7^ n — 
involves an analysis in terms of space quantization which is beyond the scope 
of the present study. 
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III. HYDROGEN MOLECULE ION 



A. Formalism 

A hydrogen molecule ion, H2, consists of two proton nuclei and one 
electron. We assume, in adiabatic approximation, the protons located at fixed 
positions z = ±c on the molecular axis (see Fig. 4). The term "Coulomb 
oscillator" denominates again the motion of a point electron, now along either 
the line through the protons {z axis) and with turning points ±C, or along 
the perpendicular line through the midpoint (y axis) with turning points ±B. 
The molecule ion's total energy E at any position on the axis, — C ^ ^ +C, 
must equal the potential energy at the turning point C, 



E = -mv^ - - r H = — \ . (7) 

2 \z + c\ \z-c\ 2c C + c C-c 2c ^ ^ 

For electron positions beyond the protons, z > c, this gives an electron speed 



, 2e / ^ C 



For positions between the protons, < 2; < c, the corresponding speed is 



, 2e / c C 



The speed expressions will be used in the action integral, 
1 /■ , 2m 





with outer contribution 



2m 

■^out — ~ 



and inner contribution 



/2m f^'^ 
p^dz = y v{z)dz = = Aout + Ain, (9) 



j Vout{z)dz + J Vout{z)dz 



(10a) 



An = I Vin{z)dz. (10b) 
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Here is a fold-out factor to be specified below. 

If the conditions are such that the electron swings along the z axis through 
the midpoint 0, that is, f j„(0) > 0, then there exists also an oscillation along 
the y axis with lateral speed m, having the same total energy, 

E = -mu^ H = H . (11) 

2 v^2^ 2c y/WT^ 2c ^ ^ 

Equal energy at the axial and lateral turning points, E{C) = E{B), Eqs. (7) 
and (11), determines the latters' geometric dependence. 



c^ 



5= + — -3c2. (12) 



Solving Eq. (11) for the lateral speed. 



_ 2e / 1 1 

^7^\l V7T^ ~ 7WT^' ^^^^ 

provides the integrand of the action integral over a lateral oscillation, 



2m 



-B 



B 



Subtraction of the protons' mutual repulsion from the total energy E of 
axial or lateral motion, Eqs. (7) and (11), gives the electronic energy. 



e2 



Eei = E- -, (15) 

in its dependence on the proton separation, R = 2c. This brackets the 
molecular problem with known atomic results, Eq. (6), in the limits of i? = cxo 
(free H atom, Z = 1) and R = (free He~^ ion, Z = 2). For those cases, as 
well as any proton-proton distance R between, we keep the action constant, 

A = + Ay = nh. (16) 

Equation (16) is the Einstein quantization conditiorP^' — a generalization of 
Sommerfeld's quantization over separable variables — where the quantum sum 
equals the sum of action integrals over topologically independent paths in 
phase space 
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Here we treat the molecule ion only in its ground state, n — 1. Analytic 
solutions of the action integrals. Eqs. (lOab) and (14), are complicated due 
to elliptic functions. We therefore integrate numerically and visualize the 
integrals by the area under the corresponding speed curves. The bold curve 
in Fig. 5 shows the axial electron speed v{z) for a far proton separation, 
R = 6rB- The electron, in its semiclassical motion, then oscillates only about 
(and through) the right proton. The area under the speed curve, Eq. (Sab), 
proportionally represents the ground state's unity of action, m ^ v{z)dz = (f)h, 
with a fold-out factor = 2 in analogy to the free-atom case, Eq. (4). The 
thin curve shows, for comparison, the axial electron speed in a free H atom — 
familiar from Fig. 3 — centered at the same proton position, +c. The pull 
from the left proton (at — c) on the oscillating electron can be seen by the 
distortion of the speed function v{z) and the redistribution of the area under 
the curve. 

When, with closer proton separation R, as in Fig. 6, the electron swings 
past the midpoint 0, then the single-cusp speed curve v{z) from Fig. 5 — 
akin in shape to letter A — becomes double-cusped (akin to letter M) and 
symmetric with respect to the bisector {y axis). Now there is also a lateral 
oscillation with speed u{y) having the same total energy E. The equality of 
E in both cases can be seen in Fig. 6 by the equahty of axial and lateral speed 
at the midpoint, ^;(0) = u{0) — a position where the electron experiences the 
same potential in either case. For convenience the lateral speed u{y), though 
perpendicular to the proton axis, is displayed in Fig. 6 together with the axial 
speed v{z). The lateral speed curve u{y), drawn dotted, is readily recognized 
by its dome shape (fl). At the bifurcation value of the proton separation, R, 
where the electron starts swinging though the midpoint 0, the axial speed 
curve v{z) changes from its one-centered A shape to a two-centered M shape. 
The area under the axial speed curve then abruptly doubles, M(i? — 6) ^ 
2A{R + d), upon a very small change in proton separation, 6 R. In order 
to keep the action integral continuous at R, the sudden area doubling is 
compensated by a corresponding doubling of the fold-out factor from (f) — 2 
for R > R to (j) = 4 for R < R. The same fold-out factor, = 4, must be 
used for the lateral action integral Ay, Eq. (14), as will become clear shortly. 

The unity of action, A = Ih, is visualized again in Fig. 6. To this end we 
compare the right half of the bold M-shape curve of axial electron speed in 
the molecule with the thin curve A (if) of the axial speed in a free if atom 
positioned at the right nucleus, +c. Due to attraction from the left nucleus. 
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the right wing of the M curve is smaller than that of A(iJ) by the area of 
lobe L. On the other hand, the left flank of the A{H) curve that extends over 
the negative z axis is smaller than the left wing of the lateral speed curve fl 
by the area of slice S. The area under the free-atom curve is then 

A{H) ^ Im + L + i n -5. (17a) 

The areas of lobe and slice are comparable, 

L^S. (17b) 

When the tiny notch to the right of the saddle point of M is taken into 
account, then the approximations (17ab) become equations and combine to 

M + n = 2A(if). (18) 

The area under both the axial and lateral speed curves is thus four times 
the area under one wing of the free-atom curve, M+fl = 4 x |A(if). Since 
the latter represents one quantum of action, h, the combined area M+fl 
visualizes its double fold-out, = 4. 

With very close proximity of the nuclei, R ^ 0, the crests of the M curve 
start merging while its saddle point, f (0), keeps rising. In the R = limit of 
fusing nuclei the axial electron speed becomes that of a free He~^ ion, M— »• 
A{He~^), familiar from Fig. 3. Concurrently, the lateral speed curve fl rises 
at its peak, u{0) = f (0), and narrows at its base until it, too, turns into the 
speed curve of the free He~^ ion, fl —>■ A{He~^). In the R = limit the three 
curves merge, M(0) = n(0) = A(iJe+). 

The results of the Coulomb-oscillator approach will be compared with 
another semiclassical calculation of by Strand and Reinhardt.!^ These 
authors, like Pauli,'^ separate the equation of motion in spheroidal coordi- 
nates, ^ = (r_|_ + r_)/2c , T] = (r+ — r_)/2c and ip by virtue of the constants 
of the motion: total energy E, angular momentum M about the z axis, and 
a component of the bifocal Runge-Lenz vector, ^Ic^ Here r+ (r_) is the 
distance of the electron from the nucleus at +c (— c). Strand and Reinhardt 
(SR) solve the ensuing one-dimensional differential equations with classical 
Poisson-bracket techniques. They find the electron's trajectories condition- 
ally periodic!^ and regionally confined due to restrictions from E, M and fic- 
However, unlike Pauli, who used Sommerfeld quantization, SR employ the 
Einstein-Brillouin-Keller (EBK) quantization conditions. 



9 




Pjdj = (uj + -)h, j = C, V 



(19ab) 



and 



n^h. 



(19c) 



The background of EBK quantization touches on the foundations of clas- 
sical and quantum mechanics.^ For the present purpose its essential rationale 
may be summarized as follows: A semiclassical treatment envolves turning 
points of radial, or other librating motion. Any tunneling through "forbid- 
den" regions of negative kinetic energy is ruled out. Viewed in terms of the 
quantum mechanical WKB (Wentzel-Kramers-Brillouin) approximation, the 
hard reflection of a wavefunction at a turning point corresponds to a so-called 
"loss" of phase (phase shift by vr) compared to the soft reflection caused by 
tunneling (phase shift by 7r/2 ). This shortcoming can be remedied with 
EBK quantization conditions by addition of a value of 1/4, for each libra- 
tional turning point, to the corresponding quantum number. Such is the case 
for the electron's elliptical and hyperbolic librations in the above quantiza- 
tion, Eq. (19ab), but not for a rotation about the z axis, Eq. (19c). Strand 
and Reinhardt call these quantization conditions "primitive" to distinguish 
them from more sophisticated ones, specified below. 

A semiclassical treatment of a free H atom with EBK quantization has 
recently been presented in these pages.l^ In this case the isotropic symmetry 
permits a separation of variables in spherical coordinates, r, 9 and y?, and 
the EBK quantization conditions are like Eqs. (19abc) except for j = r, 9. 
The atomic ground state is characterized by the quantum numbers (n^, n^, 
fT-if) = (0, 0, 0). Accordingly, the action in the atom's ground state, Ai = Ih, 
is attributed only to tunneling at the radial and latitudinal turning points 
(the former being the nucleus). Applying EBK quantization to the ground 
state of H2, denoted lsJ2g in molecular spectroscopy, SR likewise assign the 
quantum numbers {n^,nri,n^) = (0,0,0). 

B. Results 

Energies of the Coulomb oscillator (CO), in adiabatic dependence on the 
proton separation, are listed in Appendix B and shown in Fig. 7 in com- 
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parison with exact quantum mechanical (QM) results and the semiclassical 
calculation by SR.!^ The lower part of Fig. 7 shows the electronic energy 
Eei{R) of the ground state, IsEg. At large proton separations, R > 6 r^, 
both the SR calculation (circles) and the CO approach (crosses) agree ex- 
cellently with the exact QM values (curve). This is the situation where the 
electron stays near one nucleus (see Fig. 5). As Fig. 7 further shows, such 
agreement ceases once the classical electron motion leads beyond the molec- 
ular bisector, which happens for proton separations below the bifurcation 
value, R < R ^ 5.57 rs (see Fig. 6). Interestingly, the deviations of CO 
and SR from QM are opposite over the entire range. The SR results are 
discontinuous at a certain proton separation, R* ^ 1.38 r^. Remarkably, at 
(or near) that value, R*, the CO energy crosses the curve of the QM solution. 

Adding to the electronic energy Egi the proton-proton repulsion gives the 
total energy E, Eq. (15). The middle part of Fig. 7 shows by the dashed 
curve the exact total energy E{R) of the IsS^ ground state, obtained from 
QM and, by symbols, the corresponding CO and SR values. The solid dot at 
the minimum of the curve shows the QM equilibrium energy, in agreement 
with experiment, Eq = —1.20 Ry, and the equilibrium internuclear distance, 
Ro = 2.00 tb- The CO energy (x) comes out too low, due to the inaccuracy 
of its Eel, with an equilibrium value Eq{CO) ~ —1.38 Ry at Ro{CO) ~ 2.5 
tb- Conversely, the SR energy (□) comes out too high with Eq{SR) ~ —1.05 
Ry at Ro{SR) ~ 5 r^. Since the CO and SR results deviate about equal and 
oppositely from QM, their average {(}) is close to the exact values with a 
minimum of ^o[|(CO + SR)] « -1.19 Ry at i?o[|(CO + SR)] ^ 2.5 r^. 

The molecular binding energy is the difference of E in the molecule and in 
the constituting atoms, here, EqIH^)— E{H). The energy of a free hydrogen 
atom, E{H) = —1 Ry, is indicated in Fig. 7 by the fine horizontal line. 
Both semiclassical treatments, CO and (barely) SR, yield molecular binding 
energies with negative values and thus a stable molecule ion. Why are 
they more successful than the early attempts, in the 1920s, by Pauli,^ and 
independently Niessen,'^with the Sommerfeld quantization conditions of the 
old quantum theory? 

For the same reason that Sommerfeld had excluded the angular quan- 
tum number / = for the H atom — avoidance of electron collision with 
the nucleus — both Pauli and Niessen excluded electron motion in the nu- 
clear plane of the H2 molecule ion. They then found the lowest admissible 
quantum state to be (n^, n^, n^) = (0, 1, 1), denoted 2pllu in molecular spec- 
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troscopy, with Eo{P,N) = -0.52 Ry at Ro{P,N) = 5.53 ±0.01 r^, depicted 
by the triangle in the top part of Fig. 7. The QM energy^ of that quantum 
state is Eo{2pIlu) = —0.27 Ry at Rq{2pIIu) ~ 8 r^, marked by the sohd dot at 
the minimum of the dotted curve. Since both these energy values are higher 
than the ground state of a free hydrogen atom, they give rise to positive 
molecular binding energies and thus to spontaneous dissociation, if^(2pnu) 
^ H + H^. Qualitatively, Pauli's and Niessen's finding of energetic insta- 
bility is borne out by quantum mechanics for this exited state of H2 (Pauli's 
argument^ about "dynamical stability" notwithstanding). The deviation of 
their historical value (A) from the (dotted) QM curve is remarkably small — 
comparable to those of the CO and SR results for the ground state. Pauli's 
and Niessen's misfortune, though, was that they misinterpreted their result 
as the molecule ion's ground state — an assessment with fateful consequences 
in the development of quantum theory. 

C. Discussion 

Why do the semiclassical results of the Coulomb oscillator and of SR's 
quantization deviate from the QM solution of the H2 molecule ion? Strand 
and Reinhardt explain the deviation of their primitive quantization from QM 
with effective potential barriers arising from constrictions due to conserva- 
tion of both the energy E and the bifocal Runge-Lenz component flc (the 
angular momentum vanishes for the ground state, M = 0). The most drastic 
consequence of those barriers is the discontinuity of E^i at R* (see Fig. 7) 
and the large deviations at closer proton separation, R < R*. While the sim- 
ulation of quantum mechanical tunneling beyond the semiclassical turning 
points of librations is adequately achieved by EKB quantization under the 
isotropic symmetry of a free H atom^ Eq. (19) is less successful under the 
lower symmetry of H2 ■ When SR remedy the situation with "unified" semi- 
classical quantization conditions, then E^i agrees, for all practical purposes, 
with QM. Those unified quantization condition, going well beyond Eq. (19), 
are sophisticated in their dependence on Egi, Qc, and the hyperbolic turning 
points ?7_|_. They will not be discussed here. 

The reason for the deviation of the CO results from the QM values is the 
neclect of the electron's wave nature in the underlying quantization condition, 
Eq. (16). In proposing his wave hypothesis de Broglie^^ already showed 
that the quantization condition of the Bohr model. An = nh, is equivalent 
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to n standing waves along the nth Bohr orbit. If s denotes the position 
along the Bohr orbit, then the de Broglie wave can be expressed as w{s) = 
sm[27ian{s)/h], with the variable an(s) = (A„/S'„) Jq ds' along the orbit's 
circumference Sn = STrr^. A generalization gives for the Coulomb oscillator 
of a free H atom in the ground state {n = 1), its de Broglie wave as 



and the speed v{z') from Eq. (2). This de Broglie wave, shown by the line 
curve in Fig. 8 for a free H atom positioned at +c, has a node at the nucleus 
and at each turning point. By the above characterization, those turning 
points are "soft," and it is their softness that ensures the exact energy of the 
free atom. However, when Eq. (20ab) is applied to H2 for proton separations 
beyond the bifurcation value, R > R, with axial speed from Eq. (Sab) 
and integration away from the occupied nucleus, then the de Broglie 

wave is found to be "truncated" (no nodes) at both the outer and inner 
turning point (see Fig. 8, small circles). Those turning points are "hard" 
and give rise to incorrect energies. Qualitatively, an augmentation of the 
truncated de Broglie wave with (exponential) "tunneling tails," determined 
by the negative kinetic energy in the classically "forbidden" region, would 
"soften" the turning points. This would give rise to an effective far turning 
point farther out, Ce// > C, and accordingly raise the CO energy, Eq. (7), 
toward the QM result. 

When the proton separation is below bifurcation, R < R, then Eq. (20b) 
should be integrated from the midpoint in both the axial and lateral direc- 
tion rather than favoring one nucleus with the node of the de Broglie wave. 
Again, the de Broglie wave is found to be truncated at the axial and lateral 
turning points, C and B, respectively (see Fig. 9). An exception exits for 
the proton separation R*. The de Broglie wave, shown in Fig. 10, then has 
a minimum at both turning points, w{C) = w{B) = —1, according to action 
values of = |/i and Ay = ^h. Such turning points seem to be "benign" — 
reminiscent of the soft turning points in the free-atom case — and cause the 
CO energy Eei{R*) in Fig. 7 to be exact. At still smaller proton separation, 
R < R*, the de Broglie wave is truncated again (not shown). The limit R = 
corresponds to a Coulomb oscillator in the free He~^ ion which, like in the 



w{z) 



sin[27ra(^)//i] 



(20a) 



with 




(20b) 
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free H atom for R = oo, has de Broglie nodes at the turning points and an 
exact energy value. 

If the explanation that the CO energy Eei{R) deviates from the QM 
curve because of the neglect of wave effects is valid, then this sheds new light 
on the SR results. The opposite sign of the CO and SR deviations then 
suggests that SR's primitive EBK quantization, while appropriate for a free 
atom, simulates too much wave effects under the lower symmetry of the H2 
molecule ion. However, the average of both those semiclassical quantizations 
seems to be an excellent compromise, as evidenced by the close agreement of 
the corresponding total energy (0) with the (dashed) QM curve in Fig. 7. 

In conclusion, semiclassical quantization can rise to Gutzwiller's chal- 
lenge and "explain" the chemical bond in the paradigm molecule,'^ ifg^, by a 
combination of classical mechanics, quantization, and moderate wave effects. 
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By Eq. (4) the action integral of the atomic Coulomb oscillator is 




For a comparison with integral tables we change notation to x = z and use 
the abbreviation a = — l/F. Then 
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V. APPENDIX A: QUANTIZATION 




Vx 



X 



dx 



(21) 



with X = ax^ + X. Integration tables give 




(22) 
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The first term on the rhs, evaluated at the hmits of the Jp integration, 
vanishes. The last integral in Eq. (22), tabulated as 



/ 



'^'^ {-Vf) arcsin(l - 2x/r), (23) 



and evaluated at the limits, x — and x — F, contributes 

— [arcsin(l) — arcsin(— l)]v^=-7rv^. (24) 

Combining Eqs. (4'), (22) and (24), together with a fold-out factor (f) — 2, 
gives the action integral, to be equated with the Sommerfeld quantization 
condition, 

^„ = -^\/2Zm-(-7r\/f) = nh. (25) 
(p 2 

We square Eq. (25) and solve for the amplitude of the quantized Coulomb 
oscillator. 



r„ = 2^r^^ (26) 



in terms of the Bohr radius r^. 



VI. APPENDIX B: DATA 



TABLE I. Electronic ground-state energy E^i for various nuclear separations 
R of the H2 molecule from quantum-mechanical calculations (QM, Rcf 21), 
the "primitive" semiclassical quantization of Strand and Reinhardt (SR, Ref. 
15), and the present Coulomb-oscillator approach (CO). 
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K [a.u.) 


QM (Ry) 


0.00 


-4.00 


0.25 


-3.80 


0.50 


-3.47 


0.75 


-3.11 


1.00 


-2.90 


1.25 


-2.68 


1.50 


-2.50 


1.75 


-2.34 


2.00 


-2.21 


2.50 


-1.99 


3.00 


-1.82 


3.50 


-1.69 


4.00 


-1.59 


5.00 


-1.45 


6.00 


-1.36 


8.00 


-1.26 


10.00 


-1.20 


12.00 


-1.17 



SR (Ry) 


CO {Ry 


-4.00 


-4.00 


-4.75 


-3.35 


-4.27 


-3.12 


-3.76 


-2.93 


-3.34 


-2.81 


-3.07 


-2.67 


-2.19 


-2.55 


-2.09 


-2.45 


-1.99 


-2.36 


-1.79 


-2.18 


-1.66 


-2.02 


-1.56 


-1.87 


-1.49 


-1.75 


-1.40 


-1.54 


-1.33 


-1.38 


-1.25 


-1.27 


-1.20 


-1.20 


-1.17 


-1.17 
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VII. FIGURE CAPTIONS 



Fig. 1. Partial trajectory of an extranuclear orbit XX' and of a pene- 
trating orbit PP' through nucleus N . The dotted line S shows the major 
symmetry axis of XX' . 

Fig. 2. Line orbit of a Coulomb oscillator with nucleus at origin and 
turning points at ±r. 

Fig. 3. Axial electron speed v vs. position z of an electron in Coulomb 
oscillation for the ground state, n = 1, of a hydrogen atom H {Z = 1, solid 
curve) and a helium ion He'^ {Z — 2, dashed curve). The area under one 
wing of each curve represents the radial action, A—lh. 

Fig. 4. Axial Coulomb oscillation in an H2 molecule ion between ax- 
ial turning points ±C and through nuclei at ±c; perpendicular oscillation 
between lateral turning points ±B and through midpoint 0. 

Fig. 5. Axial speed v vs. position z in the ground state of an H2 molecule 
ion (bold curve) and, for comparison, of a free H atom (thin curve). Circles 
indicate the axial the positions of the nuclei, here with a large separation. 

Fig. 6. Axial speed v vs. position z in the ground state of an H2 
molecule ion (M-shaped bold curve) and, for comparison, of a free H atom 
(A-shaped thin curve, centered at the right nucleus). The Pi -shaped dotted 
curve, centered at 0, shows the lateral speed u vs. the perpendicular position 
y in the molecule. Circles indicate the axial positions of the nuclei, here with 
a small separation, R = 2 rs- 

Fig. 7. Dependence of energies of H2 on internuclear distance R. Bottom: 
Electronic energy E^i from quantum mechanics (QM, solid curve), primitive 
semiclassical quantization by Strand and Reinhardt (SR, Q) ^^'^ ^^e present 
Coulomb-oscillator approach (CO, -|-). Middle: Total energy E of the IsEg 
ground state from QM (dashed curve) with minimum (•), values by SR (□) 
and CO (x) and their average (o), energy of a free H atom (dashed horizontal 
line). Top: Total energy E of the 2pllu state from QM (dotted curve) with 
minimum (•), and historical value by Pauli and Niessen (A). 

Fig. 8. de Broghe wave of the free H atom Coulomb oscillator (curve) and 
of the Coulomb oscillator (small circles) for the same proton separation 
as in Fig. 5, R = 6 tb- Large circles indicate the axial positions of the nuclei. 

Fig. 9. de Broglie wave of the H2 Coulomb oscillator along axial distance 
OC from Fig. 4 (right side of graph) and lateral distance OB (left side). 
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Large circles indicate the axial positions of the nuclei, here with the same 
separation as in Fig. 6, i? = 2 r^. 

Fig. 10. Same as Fig. 9 but for the proton separation R* — 1.38 
where the semiclasssical energy is exact. 
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